Effects of periodic potentials on the critical velocity 
of superfluid Fermi gases in the BCS-BEC crossover 



Gcntaro Watanabc, 1,2 Franco Dalfovo, 3 Lev P. Pitacvskh, 3,4 and Sandro Stringari 3 

'Asia Pacific Center for Theoretical Physics (APCTP), POSTECH, 
San 31, Hyoja-dong, Nam-gu, Pohang, Gyeongbuk 790-784, Korea 
2 Nishma Center, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan 
3 INO-CNR EEC Center and Department of Physics, University of Trento, 38123 Povo, Raly 
Jt Kapitza Institute for Physical Problems, 119334 Moscow, Russia 
(Dated: July 26, 2011) 

We study the effects of an external periodic potential on the critical velocity of a superfluid 
Fermi gas in the crossover between the Bardeen-Cooper-Schrieffer (BCS) phase and Bose-Einstein 
condensation (BEC). We numerically solve the Bogoliubov-de Gennes equations to model a three- 
dimensional (3D) gas of ultracold atoms in the superfluid phase flowing through a ID optical lattice. 
We find that when the recoil energy is comparable to the Fermi energy, the presence of the periodic 
potential reduces the effect of pair-breaking excitations. This behavior is a consequence of the 
peculiar band structure of the quasiparticle energy spectrum in the lattice. When the lattice height 
is much larger than the Fermi energy, the periodic potential makes pairs of atoms to be strongly 
bound even in the BCS regime and pair-breaking excitations are further suppressed. We have also 
found that when the recoil energy is comparable to or larger than the Fermi energy, the critical 
velocity due to long-wavelength phonon excitations shows a non-monotonic behavior along the 
BCS-BEC crossover. 

PACS numbers: 03.75.Ss, 03.75.Lm, 67.85.De 



I. INTRODUCTION 

Ultracold atom gases in optical lattices have been con- 
tinuously attracting great interest for the last ten years 
[H-Q . Recent developments in the field of ultracold atom 
gases provide a new research arena in the physics of 
quantum fluids: by using Feshbach resonances of ul- 
tracold Fermi atoms, one can study the crossover from 
the Bardeen-Cooper-Schrieffer (BCS) phase to a Bosc- 
Einstein condensate (BEC) of molecules In the cur- 
rent research frontier, superfluidity of ultracold Fermi 
gases in optical lattices is very intriguing problem, which 
has interesting connections with similar issues in solid 
state physics, nuclear physics, and astrophysics. 

The critical velocity of supcrfiow due to energetic in- 
stability is one of the most important properties of su- 
pernuids, which has been pioneered by Landau [f|. If 
the velocity of superfiow exceeds some critical value, the 
kinetic energy of the superfluid can be dissipated by 
creating excitations 043- ^ n uniform superfluid Fermi 
gases in the BCS-BEC crossover, excitations which cause 
the energetic instability are of two types: fermionic 
pair-breaking excitations in the BCS regime and long- 
wavelength phonon excitations in the BEC regime [l,[lojj. 
In the unitary regime both mechanisms are suppressed 
and the critical velocity shows a maximum value 0-[ll|. 

Recently, effects of periodic potentials on the critical 
velocity of Fermi supcrfluids has been studied experimen- 
tally [12| ■ This experiment has stimulated theoretical in- 
vestigations of this problem [l3l - [Tfl |. Most of them has 
focused on the BCS regime in tight-binding approxima- 
tion (l3l - [l5l | . The purpose of the present work is to obtain 
an understanding of the critical velocity from a unified 



point of view covering all regions along the BCS-BEC 
crossover and both the strong and weak lattice regime. 
To this purpose, we use the Bogoliubov-de Gennes (BdG) 
equations. This theory accounts for both types of exci- 
tations which are relevant in this problem. In our previ- 
ous work [l6| we already used it for a gas at unitarity; 
here we extend the calculations in order to explore the 
whole crossover region. As a main result, we find that, 
when the lattice height is comparable to or much larger 
than the Fermi energy, the periodic potential reduces the 
effect of pair-breaking excitations. This is due to the 
periodic structure of the quasiparticle energy spectrum 
in the Brillouin zone and the formation of the bound 
molecules induced by the lattice. Another main result is 
that when the recoil energy is comparable to or larger 
than the Fermi energy, the critical velocity due to long- 
wavelength phonon excitations shows a non-monotonic 
behavior along the BCS-BEC crossover. These effects 
are unique for Fermi supcrfluids in periodic potentials 
and do not exist in the case of single barrier potentials 

EMI- 

This paper is organized as follows. In Sec. [Til we ex ~ 
plain the basic formalism employed in the present work. 
Then we show the results in Sec. IIIII Finally, summary 
and outlook are given in Sec. IIVI 



II. BASIC FORMALISM 

We want to study the effect of the periodic potential 
on the Landau critical velocity of Fermi supcrfluids in the 
whole BCS-BEC crossover, in situations where the Fermi 
energy is larger or smaller than the lattice height. For 
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this aim, we need to use a theoretical framework which 
can account for the formation of bound molecules induced 
by the periodic potential, which is important when the 
lattice height is larger than the Fermi energy [20| ; 
the same formalism must also account for pair-tunneling 
processes, which are important on the BEC side of the 
resonance [2(| Hl| . A suitable approach consists of the 
numerical solution of the Bogoliubov-de Gennes (BdG) 
equations [22j]: 
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where Ui and Vi are quasiparticlc amplitudes and the 
corresponding cigen-encrgics. The single-particle hamil- 
tonian is H'(r) = — fi 2 V 2 /2m + V cx t — where m is the 
atom mass and T4 x t (r) is the external potential. The or- 
der parameter (or gap parameter) A(r) and the chemical 
potential ^, appearing in Eq. ([I}, are variational param- 
eters determined from the gap equation, 



A ( r ) = -.9 5>*( r K( r ) 
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together with the constraint 
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enforcing the conservation of the average density no. 
Here g is the coupling constant for the contact in- 
teraction and V is the volume of the system. The 
BdG eigenfunctions obey the normalization condition 
/ d 3 r [u*(r)uj(r) + v*(r)vj(r)] = Sij. Finally, the energy 
density e can be calculated as 

e =V /rfr^^A-^I^Wr + A^rKW^W]. (4) 

In the present study, we consider a three-dimensional 
supcrfluid Fermi gas, which is uniform in the x and y di- 
rections and subject to a one-dimensional optical lattice 
along z: 

V ext (z) = sE R sin 2 q B z = V sin 2 q B z . (5) 

Here Vb = sE R is the lattice height, s is the laser inten- 
sity in dimensionless units, E R = h 2 q B /2m is the recoil 
energy, q B = n/d is the Bragg wave vector, and d is the 
lattice constant. For practical reasons, throughout this 
paper, we set s = 1 except for special cases, which we 
shall mention explicitly. The ratio between the Fermi 
energy and the lattice height is then varied by changing 
the average density of the gas. 

In the presence of a supercurrent with wave vector Q = 
P/h moving in the direction of the periodic potential, one 
can write the gap parameter in the form 



where A(z) is a complex function with period d. There- 
fore, from the gap equation, we see that the eigenfunc- 
tions of Eq. ([T|) must have the Bloch form Uj(r) = 
Ui(z)e l Q z e l1i ' T andi>i(r) = Vi(z)e~' l< ^ z e lk ' r . The wave vec- 
tor k z lies in the first Brillouin zone (i.e., — q B < k z < q B ) 
and Ui and Vi are periodic in z with period d. We also 
define the quasi-momentum P c d gc and quasi-wavenumbcr 
Qedge at the edge of the Brillouin zone for P and Q as 
-Pcdgc = ftQedge = %b/2. (Note that the edge of the first 
Brillouin zone for P is at P c d gc = <7b/2 while that for k z 
is at q B .) This Bloch decomposition transforms Eq. ([T]) 
into the following BdG equations for Ui and vf. 
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where 



H' Q (z) 
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(8) 

Here, the label i represents the wave vector k as well 
as the band index. In order to remove the ultravio- 
let divergences in the BdG equations with contact po- 
tentials, we use the regularization scheme proposed by 
Refs. [23|, Since we need to calculate the second 

derivatives of the energy with respect to the density and 
the quasi-momentum, we use large values of the cutoff 
energy Ec, especially in the BEC side, where the size of 
the pair is much smaller than the average inter-atomic 
distance (details are described in Sec. IIII[) . 

As discussed in Refs. 0, [lfj, the energetic instabil- 
ity of superfluids of dilute Fermi gases can be caused by 
two processes [25| : the creation of long- wavelength super- 
fluid phonon excitations or fermionic pair-breaking exci- 
tations. The critical velocity by the former process can 
be determined by the hydrodynamic analysis of the exci- 
tations 0, [H, [28l - [30| . Starting from the continuity equa- 
tion and the Euler equation, and linearizing with respect 
to the perturbations of the density and the velocity fields, 
we obtain the dispersion relation of the long-wavelength 
phonon, 
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Here, Tiuj and q are the energy and the wavenumber of 
the excitations, rio and P are the average density and 
the quasimomcntum of the superfluids. The energetic 
instability occurs when u(q) becomes negative: 
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A(r) = e l2Qz A(z), 



(6) 



In practice, we calculate the energy density e(no,P) for 
a given uq and P from Eq. (jl]) using the solution of the 
BdG equations (JT]). Then the critical quasimomentum P c 
at which the energetic instability occurs is determined by 
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and k F = (37r 2 no) 1 ^ 3 are the Fermi energy and mo- 
mentum, respectively, of a uniform noninteracting Fermi 
gas of density no. For each value of Ep/En, we solve 
the BdG equations for several values of the parameter 
l/k F a s along the crossover from the BCS to the BEC 
side, namely l/k F a s = — 1, —0.5, 0, 0.5, and 1, where a s 
is the s-wave scattering length of atoms. 

In the x and y directions, we assume periodic boundary 
conditions with period L± = 127rfcp 1 . We set the cutoff 
energy E c as follows: for E F /E R = 2.5, E c = 40i?F 
in the BCS side and Ec = IOOEf at unitarity and in 
the BEC side; for E F /E R = 1, E c = 50£ F in the BCS 
side and Ec = 100-Ef at unitarity and in the BEC side; 
for E F /E R = 0.1, E c = 250£ F in the BCS side, E c = 
ZbOE F at unitarity, and E c = 500-E F in the BEC side. 



FIG. 1: (Color online) Density profile n(z) at P — in the 
BCS-BEC crossover for Ef/Eb, — 1 and s = 1. Here, n(z) is 
normalized by the average density no 



Eq. (flO|) [3lj . We finally obtain the critical velocity v c 
from 

1 ( de\ 

Vc = — 7p • 11 

n \dPJ Pc 

On the other hand, the critical velocity due to the 
pair-breaking fermionic excitations can be determined by 
looking at the quasiparticle energy spectrum e^. The en- 
ergetic instability by the pair-breaking excitations occurs 
when some quasiparticle energy starts to be negative: 

e t < . (12) 

From Eq. (|11[) evaluated at the critical quasi-momentum 
determined by this condition, we obtain a critical veloc- 
ity for the pair-breaking excitations. The actual critical 
velocity of the system is the lowest between the ones ob- 
tained from the above two conditions [32[. We finally 
note that the gas becomes unstable also when some ex- 
citation energy starts to have a non-zero imaginary part. 
This corresponds to a dynamical instability, which causes 
an exponential growth of the amplitude of the perturba- 
tion. To address the problem of dynamical instability, 
short- wavelength bosonic excitations should be also prop- 
erly included. This is beyond the scopes of the present 
work, in which we instead focus on the energetic insta- 
bility. Results of the critical velocity for dynamical in- 
stability due to long-wavelength excitations are given in 
Appendix O 



III. RESULTS 

We study the three cases of E F /E R = 2.5, 1, and 
0.1 with a fixed value of s = 1 except for a few cases 
which we shall mention explicitly. Here E F = h 2 k F /(2m) 



A. Density profiles and gap parameter 

In Fig. [TJ we show the density profile n(z) = 

2J2i \vi( z )\ of Fermi atoms at P = along the BCS- 
BEC crossover for E F /E R = 1. Moving from the BCS 
regime (l/k F a s = —1) to the BEC regime (l/k F a s = 1), 
the density n(z) becomes more inhomogeneous. This be- 
havior is consistent with the fact that the compressibility 
of a uniform Fermi gas is known to increase monotoni- 
cally with l/k F a s . We observe the same qualitative be- 
havior for the other values of E F /E R . 

In Fig. [21 we show the amplitude |A(z)| of the order 
parameter at P — for different values of l/k F a s and 
E F /E R . From this figure one can see that, especially 
in the BCS regime (l/k F a s = —1), the order param- 
eter |A| is enhanced when the Fermi energy is smaller 
than the lattice strength, as shown by the blue curves for 
E F /E R = 0.1 (i.e., E F /V = 0.1). We understand this 
fact as due to the formation of bosonic molecules induced 
by the external periodic potential. This process is indeed 
expected to become significant when the lattice is strong 



B. Critical velocity 

Before presenting the numerical results for the criti- 
cal velocity, let us discuss the conditions upon which the 
periodic potential can produce significant effects on the 
behavior of long-wavelength phonons and pair-breaking 
excitations. For long-wavelength phonons the condition 
is E F /E-R > 1 or ,s ^> 1. In fact, in the opposite case, 
E F /E K <C 1 and s < 1, the right-hand side of Eq. (TlU)) , 
which coincides with the sound speed at P — 0, is much 
smaller than the maximum value of the left-hand side 
of Eq. (flOl). which is of order q^jm (see also Fig. 5 in 
Ref. [l||). Thus Eq. (TO]) is satisfied at P ~ and 
the critical velocity is basically determined by the sound 
speed at P = 0, which is also close to the sound speed 
in the uniform system. Consequently, the critical ve- 
locity due to long-wavelength phonons is almost unaf- 



4 



r fo 



0.5 
0.4 
0.3 
0.2 
0.1 


1.2 



(cl) l/k F a s = -1 



r ^ 


1 


&3 




\ 


0.8 




0.6 


< 






0.4 




0.2 








r ^ 



< 




2 
1.5 
1 

0.5 




FIG. 2: (Color online) Amplitude |A(z)| of the gap param- 
eter at P = in the BCS-BEC crossover: l/fcpa s = —1 (a), 
(b), and 1 (c). The horizontal dotted lines show the ampli- 
tude of the gap parameter for the uniform system at the same 
value of l/fcpa s . 



fected by the presence of the lattice if Ep/En <C 1 and 
s < 1, even though the lattice height Vo is large com- 
pared to the Fermi energy Ep. On the other hand, the 
sufficient condition for pair-breaking excitations to be af- 
fected by the periodic potential due to the formation of 
bound molecules is Ep/Vo <C 1. This condition can be 
satisfied either by decreasing Ep/E^ or by increasing s. 
When Ey/En ~ 1, the critical velocity for pair-breaking 
excitations is also affected by a peculiar band structure 
of the quasi-particle spectrum. 

Our results for the critical velocity v c are shown in 
Fig. [3] for s = 1. Let us first concentrate on the results at 



high density, E-p /En = 2.5, in panel (a). The open circles 
correspond to the critical velocity for long-wavelength 
phonons, which exhibits a non-monotonic behavior. In 
particular, in the BCS regime (negative l/kpa s ) this crit- 
ical velocity is strongly reduced compared to the one in 
a uniform gas (see, e.g., Fig. 8 in Ref. This is a 

peculiar effect of the lattice. However, for the parame- 
ters of Fig. EJa), the actual critical velocity in the BCS 
regime is still given by fermionic pair-breaking excita- 
tions (filled squares). The latter are almost unaffected 
by the lattice and therefore, the overall behavior of the 
critical velocity in the crossover is qualitatively similar to 
that of a uniform gas, already discussed in Ref. [lfj: in 
the BCS regime, v c increases when approaching unitarity 
(l/kpa s = 0), because the intra-pair attraction becomes 
stronger and thus the amplitude of the gap parameter 
increases; in the opposite BEC regime (\/kpa s > 0), the 
critical velocity is given by long- wavelength phonons and 
an increase of the inter-pair repulsion leads to a larger 
sound speed and, again, a critical velocity v c increases 
towards unitarity. As a consequence, v c takes a maxi- 
mum value at l/kpa s ~ 0. 

When the recoil energy is comparable to the Fermi 
energy, the periodic potential causes qualitative changes 
in the results of the critical velocity. For Ep/Epi = 1 
[Fig-EIb)], we observe that, at l/kpa s = —0.5, the crit- 
ical velocity is given by long-wavelength phonon exci- 
tations rather than pair-breaking excitations even in the 
BCS regime. We understand this effect as mainly due to a 
peculiar band structure of the quasiparticle energy spec- 
trum. In Fig. 21 we show the lowest band of the quasipar- 
ticle energy spectrum q for the first radial branch with 
k]_ = k 2 x + k 2 y = at P = and l/k F a = -0.5. In gen- 
eral, the quasiparticle spectrum near the center of the 
Brillouin zone, at \k z \ ~ 0, is only weakly affected by the 
periodicity of the system and hence the change of with 
increasing P is close to that in the uniform system, given 
by the Doppler term Phk z /m. On the other hand, close 
to the zone edge, at k z ~ ige, the change of q with in- 
creasing P is much smaller than Phk z /m because of the 
periodicity of the Brillouin zone (e, at k z = ±<7b must 
be identical). In the case of Ep/En = 1, the minimum 
of e, is indeed located close to the edge of the Brillouin 
zone, unlike the other two cases of Ep/Eji = 2.5 and 0.1. 
Therefore, the reduction of the minimum value of ti with 
increasing P is relatively small for Ep/E^ = 1 and this 
is why we find that v c is determined by phononic instead 
of fermionic excitations in this case. 

For smaller density (Ep/E-R = 0.1), we observe a sig- 
nificant increase of v c in the whole crossover [Fig. Etc)]. 
In our previous article [l6j . we already showed that, in a 
unitary Fermi supcrfluid with Ep / Epi <C 1, the phononic 
critical velocity is almost unaffected by the lattice, if 
s < 1, and remains close to the speed of sound of a 
uniform gas with the same density 34]. This can be 
understood by recalling that phonons always have wave- 
length larger than the healing length of the supcrfluid, 
which is the order of fcp 1 or greater. When Ep/E-R <C 1, 
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FIG. 3: (Color online) Critical velocity v c of the energetic instability for Ef/Er = 2.5 (a), 1 (b), and 0.1 (c) with s = 1 in the 
BCS-BEC crossover. Open circles and filled squares show the critical velocity due to long-wavelength phonons and fermionic 
pair-breaking excitations, respectively. The horizontal dotted line in panel (c) represents the value of the sound velocity ci ' of 

uniform system at unitarity, ci°'/uF = (1 + /3) 1//2 /y / 3 ~ 0.443. The red solid lines and the black dashed lines are guides to the 
eye. 




0.2 - 

f ■ 1 ■ 1 ■ 1 ■ 1 ■ 

0.2 0.4 0.6 0.8 1 

k z /q B 



FIG. 4: (Color online) Lowest band of the quasiparticle en- 
ergy spectrum ej for P = and l/fea = —0.5. Here, we show 
the first radial branch with k\ = k^+ky = 0. For Ef/E-r = 1, 
the minimum of ej is located close to the Brillouin zone edge 
k z = qs ■ 



0.1 




10 1 1 

-1 -0.5 0.5 1 

1/ k F a s 

FIG. 5: (Color online) Critical current j c = nov c for the 
same cases in Fig. [3] Here we show the lowest value of j c , given 
either by long-wavelength phonon excitations or by fermionic 
pair-breaking excitations. The curves connecting symbols are 
guides to the eye. 



the healing length becomes much larger than the lattice 
spacing d = ir/qn, and this makes the phonons insensitive 



to the lattice itself. In the present work, we find that the 
same is true even away from unitarity, leading to a larger 
critical velocity in the whole crossover [empty circles in 
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FIG. 6: (Color online) Lowest band of the quasiparticle en- 
ergy spectrum ti for large lattice height with a = 5 and 
E F /E R = 0.1 (i.e., E F /V = 0.02) in the BCS regime at 
l/fcpa a — — 1. Here, we show the first radial branch with 
k'^ = k'i + ky — 0, which always gives the smallest values of 
e; in this case. The inset shows the amplitude |A(z)| of the 
order parameter at P = 0. The horizontal dotted line shows 
the amplitude of the order parameter for the uniform system 
at the same values of l/k F a a — — 1. 



Fig-EHc)]. The critical velocity due to pair-breaking exci- 
tations (filled squares) is also increased because a lattice 
strength Vq much larger than the Fermi energy gives a 
stronger attraction between paired atoms. 

In Fig. [5] we show the critical current j c = uqv c for the 
same cases of Fig. [3] Due to the low density at E-p /En = 
0.1, the critical current is much smaller than the other 
cases even though v c in units of vb = <Zb/w for Ep/En = 
0.1 is comparable to that of Ep/En = 1. 



C. Dependence on the lattice height 

All results shown in Figs. [T][5] have been obtained by 
fixing the lattice height s = 1 and vary the density in 
order to change the key parameter E-p /En- If we increase 
s keeping the average density fixed, the superfluid flow is 
suppressed and v c is also reduced in general. 

A systematic analysis as a function of s, which would 
be natural from the experimental viewpoint, is computa- 
tionally very demanding and is beyond the scopes of this 
work. The choice of s = 1 is not accidental, however. It 



turns out, in fact, that around s = 1 the effects of the 
lattice on the critical velocity are the most pronounced 
as far as the interplay between pair-breaking and long- 
wavelength bosonic excitations is concerned. For lower 
values of s these two types of excitations behave quali- 
tatively the same as in a uniform superfluid, as function 
of 1/kpa, being scarcely affected by the lattice, at least 
within the range of E-p /En considered in this work. On 
the other hand, at larger s the pair-breaking instability 
is quickly suppressed and long- wavelength excitations be- 
comes dominant along the crossover. 

The reason why a strong lattice prevents the pair- 
breaking processes can be understood by looking at 
Fig. [51 where we show the quasiparticle energy spectra 
at various values of P in the case of s = 5, Ep/En = 0.1 
(i.e., Ep/Vo = 0.02), and l/kpa s = —1. The spectrum 
for P = shows a quadratic dependence of k z with a pos- 
itive curvature around k z = and there are no minima at 
k z =^= 0. Even though the figure represents a case in the 
deep BCS regime, the structure of e; is consistent with 
the formation of bound pairs. In the inset of the same 
figure, we show the amplitude |A(z)| of the gap parame- 
ter at P = 0. First, we note that the minimum value of 
|A(z)| at z/d = ±1 is smaller than, but still comparable 
to the value of |A| in the uniform case, suggesting that 
the system is indeed in the superfluid phase. More im- 
portantly, one sees a large enhancement of |A(z)|, near 
z = 0, compared to the uniform system, which shows the 
formation of bosonic bound molecules. A consequence 
of this lattice induced molecular formation is that the 
energetic instability due to pair-breaking excitations is 
suppressed and does not occur at any values of P. 

For the same parameters of Fig. the energetic in- 
stability due to long-wavelength phonons instead occurs 
at P = 0.226?m7b = 0.452P c d go and the corresponding 
critical velocity at l/kpa s = — 1 is v c = 0.0662wf. We 
also find v c = 0.0429i>f at unitarity and v c = 0.0302up 
at \/kpa s = 1. This means that, in the whole crossover, 
the critical velocity at s = 5 is largely reduced compared 
to the red line in Fig. [3Jc) for s = 1 and does not exhibit 
a maximum anymore. The behavior of v c as a function 
of s at unitarity has been already discussed in Ref. [l||. 



IV. SUMMARY AND OUTLOOK 

We have studied the effects of a periodic potential on 
the Landau critical velocity of a Fermi superfluid in the 
BCS-BEC crossover. We have considered a 3D superfluid 
Fermi gas flowing in a ID periodic potential produced by 
an optical lattice. Using the Bogoliubov-de Gennes equa- 
tions, we have obtained a unifying picture both for weak 
and strong lattices and in the whole BCS-BEC crossover. 
We have found that, when the recoil energy is comparable 
to the Fermi energy, energetic instability due to fermionic 
pair-breaking excitations can be less effective as a con- 
sequence of the periodic structure of the quasiparticle 
energy spectrum. When the lattice height is much larger 
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than the Fermi energy, pair-breaking excitations are pre- 
vented because the lattice potential gives a stronger at- 
traction between paired atoms, eventually forming bound 
bosonic molecules. We have also found that, when the 
recoil energy is comparable to or larger than the Fermi 
energy, the critical velocity due to the long wavelength 
phonon excitations is drastically reduced by the lattice 
in the BCS regime leading to its non-monotonic behavior 
along the BCS-BEC crossover. 

A further interesting issue regards the possible ex- 
istence of roton-like minima in the bosonic dispersion 
curve. This excitations are obtained at low filling frac- 
tions and within a tight-binding attractive Hubbard 
model [ill, EH HH [I?} . The roton-like minima arise from 
strong charge-density-wave fluctuations. These fluctua- 
tions are expected to be less favored in our system, where 
the gas is uniform in the transverse directions. However, 
if such roton-like excitations exist also in the our case (3D 
gas in a ID lattice), they would lower the critical veloc- 
ity in the BCS regime and for strong lattices. To address 
this issue, one should use, for instance, a quasiparticle 
random phase approximation (QRPA) on top of the sta- 
tionary solution of the Bogoliubov-de Gennes equations. 
This is an interesting challenge for future investigations. 

Finally, we would like to discuss a similarity between 
the present system and nuclear "pasta" phases [36l438| 
in crusts of neutron stars. The pasta nuclei are those 
of exotic shapes such as rod-like and slab-like structures. 
In neutron star crusts, the pasta nuclei are immersed 
in background electrons and a gas of dripped neutrons, 
which is regarded to be in the supcrfluid phase. The 
setup considered in the present work resembles supcr- 
fluid neutrons in the pasta phase with slab-like nuclei, 
which are in the normal phase and provide a ID periodic 
potential for superfluid neutrons [391 ] . 
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Appendix A: Dynamical instability within 
hydrodynamic analysis 

In this Appendix, we show the results of the crit- 
ical velocity v c for dynamical instability due to long- 
wavelength excitations. This instability occurs when a 
non-zero imaginary part appears in the excitation energy 
©. The condition is 



d 2 e d 2 e 
driQ dP 2 



< 



(Al) 



Even though the long-wavelength excitations in general 
do not give the lowest critical velocity of the dynamical 
instability, this condition becomes useful at Vq ^> Ep, 
when the P-dependence of e is almost sinusoidal and 
the critical velocity dependents rather weakly on the 
wavelength of the excitations. In this case, the above 
condition gives an onset of the dynamical instability 
at P = P dgc/2, which coincides with the condition 
e(P + hq) - e(P) = e(P) - e(P - hq) for V > |A| 
This corresponds to the energy and momentum conser- 
vation for two particles decaying into two different Bloch 
states with e(P ± hq). 

In Fig. [JJ we show the critical velocity for the dynam- 
ical instability determined by the condition (|A1[) . The 
most striking feature is the large values of v c in the case 
of s = 1 and Pf/Pr =0.1. As in the case of the energetic 
instability, this is due to the fact that phonons are insen- 
sitive to the lattice when s < 1 and E-p /En <C 1, i.e., 
the healing length is much larger than the lattice spac- 
ing. Note that, even for the same value of Pf/Pr = 0.1, 
the critical velocity can be rather small provided s> 1. 
This tendency is confirmed by the results for s = 5 and 
Ef/Er =0.1 shown by crosses in Fig.[JJ 
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